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Testing non-standard CP violation in neutrino propagation 

Walter Winter a * 
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Abstract 

Non-standard physics which can be described by effective four fermion interactions may be an additional source 
of CP violation in the neutrino propagation. We discuss the detectability of such a CP violation at a neutrino 
factory. We assume the current baseline setup of the international design study of a neutrino factory (IDS-NF) for 
the simulation. We find that the CP violation from certain non-standard interactions is, in principle, detectable 
significantly below their current bounds - even if there is no CP violation in the standard oscillation framework. 
Therefore, a new physics effect might be mis-interpreted as the canonical Dirac CP violation, and a possibly even 
more exciting effect might be missed. 
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1. Introduction 

Physics beyond the Standard Model may intro- 
duce non-standard interactions (NSI) [1-5] sup- 
pressed by a higher energy scale. In general, such 
new physics is usually described by effective di- 
mension six [6-8] and eight [9,10] operators. One 
can describe the effective dimension d Lagrangian 
as a function of the non-standard physics scale A 
as 

^ = ^' W 

where A is a dimensionless coupling constant 
and O d is a dimension d operator. Thus, 
the non-standard physics will be suppressed by 
(Eewsb/A) 6 * -4 with respect to the weak interac- 
tions, where -Eewsb is the electroweak symmetry 
breaking scale. 

In this study, we focus on non-standard propa- 
gation effects in standard oscillations (SO). These 
can be phenomenologically described by neutral 
current-type NSI of the form 



in general, e™^ are complex numbers for a ^ /3, 



c a/3 

and real numbers for a 
C^l^lexp(i^)- Thus,-, 



where we define 
e M > e ™r> e ™ are pos- 
sible sources of non-standard CP violation (nsi- 
cpv) H 

They enter the propagation Hamiltonian in fla- 
vor base proportional to the matter potential 
a cc = 2^EGpN e (with N e the electron den- 
sity) in the off-diagonal elements. Since |e™| 
is very well constrained, we focus on \e™ T \ and 
|e™|, for which the current bounds are 0(0. 1) 
and 0(1), respectively [10, 14, 15]. Therefore, 
the phases (j)™ T and 0™ might be accessible by 
future experiments for large enough \e™ T \ and 
|e™|. The necessary conditions for an underly- 
ing model producing such large NSI are discussed 
elsewhere [13,16]. Since these interactions will be 
suppressed by at least a factor of A 2 (c/., Eq. (JTJ), 
it might be plausible to look for NSI-CPV in the 
best discussed neutrino oscillation experiments 
which are sensitive to the highest A-scales, such 



G 



£ NS i = -^ef a {v p YLv a ){f lp f) + h. 



(2) 



with L = 1 — 7 D , which affect the neutrino prop- 
agation in matter for / <= {e,u, d}. Note that, 
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2 In Refs. [11, 12] such NSI-CPV was discussed in the con- 
text of source and detection NSI, whereas we focus on the 
propagation effects. Note that, depending on the model, 
source and propagation NSI could be related. However, 
the simplest allowed models to induce e™ T or e™ involve 
two mediator fields (and some cancellation conditions), 
and propagation NSI are not related to source and detec- 
tion NSI [13]. 
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as neutrino factories [17,18]. 

The measurement of NSIs in neutrino factories 
has been discussed in Refs. [11,19-25]. As illus- 
trated in Refs. [24,25], the disappearance chan- 
nels and second "magic" baseline [26] are manda- 
tory for excellent NSI sensitivities. In particu- 
lar, e™ T is best measured with the disappearance 
channel, whereas e™ is best measured with the 
appearance channel. Therefore, we expect that 
the measurement of </>™ will be qualitatively sim- 
ilar to that of 6cp, whereas that of <f>™ T will have 
completely new characteristics. We use in this 
study the baseline setup of the international de- 
sign study for a neutrino factory [27], which in- 
cludes two baselines, as well as the disappearance 
channels by the measurement of the wrong-sign 
muons. 

A focus of this letter is to demonstrate that 
the discovery of NSI-CPV should be quantified 
with performance indicators similar to SO-CPV. 
In addition, the full (relevant) parameter space 
using a full simulation is discussed. Since there is 
no model-independent connection between source 
or detection and propagation NSI [13] and there 
is not yet any near detector specification for the 
neutrino factory within the international design 
study, we do not discuss source and propagation 
NSI. 

2. Method and performance indicator 

Our simulations use the GLoBES software [28, 
29] with the current best-fit values and solar os- 
cillation parameter uncertainties from Ref. [15], 
as well as a 2% error on the (constant) matter 
density profile, i.e., we expect the matter density 
profile to be known with that precision^ For the 
sake of simplicity, we use a normal simulated mass 
hierarchy. The experimental scenario we consider 
is the IDS-NF 1.0 setup from Ref. [27], which is the 
current standard setup for the "International de- 
sign study of the neutrino factory" (IDS-NF). This 
setup has been optimized within Refs. [30,31] for 
the measurement of sin 2 29i% , the neutrino mass 
hierarchy, and leptonic CP violation in the case 

3 In fact, we have checked that the impact of a larger mat- 
ter density uncertainty on the measurement is very 
small, at the level of a few percent correction. 



of standard oscillations. In short, it uses two 
baselines at about 4 000 km and 7 500 km with 
two (identical) magnetized iron neutrino detec- 
tors (MIND) with a fiducial mass of 50 kt each. 
For each baseline, a total of 2.5 • 10 21 useful muon 
decays plus 2.5 • 10 21 useful anti-muon decays in 
the straight of the corresponding storage ring is 
used, which could be achieved by ten years of op- 
eration with 2.5 • 10 20 useful muon decays per 
baseline, year, and polarity. The muon energy 
is assumed to be 25 GcV, which is sufficient 
for a detector with a low enough detection thresh- 
old [30]. The detector and systematics specifica- 
tions can be found in Refs. [27]. Note that there 
is not yet any near detector specification. We 
do not simulate the near detector explicitely, be- 
cause we do not discuss non-standard production 
or detection effects such as in Ref. [32]. As a 
small modification of the IDS-NF baseline setup, 
we do not include the emulsion cloud chamber 
for v T detection at the short baseline, since it has 
been demonstrated in Ref. [25] that it hardly con- 
tributes to the SO and NSI sensitivities if two 
baselines are used. We have checked that this 
also applies for a (hypothetical) — > v T oscilla- 
tion channel for the effects discussed in this study 
(which might be different for NSIs in the produc- 
tion process, see Ref. [12]). 

We define the sensitivity to NSI-CPV, in the 
same way as the sensitivity to SO-CPV, as the 
A% 2 with which any CP conserving solution can 
be excluded. That is, we simulate a true e™^ = 
kSal ex P(*0Q/3)> where <j>™a is the CP violating 
phase ^ {0, 7r}. In addition, we have a set of sim- 
ulated values for the SO parameters. Then we 
compute the Ax 2 for </>™^ (fit) fixed to and ir 
(CP conservation) and choose the minimum be- 
tween these two values. All the fit SO parameters 
and |e™^| (fit) are marginalized over. For the sake 
of simplicity, we do not take into account the mass 
hierarchy degeneracy. 

3. Discovery of non-standard CP violation 

In the context of NSI-CPV, the most impor- 
tant question might be for which region of the 
parameter space NSI-CPV will be discovered at a 
neutrino factory. Since the solar and atmospheric 
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NSI-CPV discovery (e£. ) NSI-CPV discovery (e£ ) 




Figure 2. Fraction of <$>™ (left) and <f>™ T (right) for which non-standard CP violation will be discovered (at 
— 9) as a function of |e™ r | (left) and |e™ | (right). The thick curves are computed for sin 2 28 13 = 0. 
The shaded regions illustrate the dependence on the (true) sin 2 26*i3 and 5cp ■ 
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Figure 1. Sensitivity to NSI-CPV in e™ as a func- 
tion of the true values of |e™ | and c/>™ (for the true 
sin 2 26*13 = 0). The different contours correspond 
to Ax 2 = 1 (dashed), 4 (thin solid), and 9 (thick). 



oscillation parameters are very well known, the 
performance will mainly depend on the true val- 



ues of |e^|, 0^g, sin 2 26»i3, and <5 C p, where </>™p 
describes the CP violation of interest. Because 
the absolute value of e™ 1 ^ suppresses the phase 
measurement, the simulated \e™p\ and <f>™p will be 
the most important parameters for the parameter 
space test, similar to sin 2 26*13 and Sep for the SO- 
CPV. We illustrate this dependence in Figure Q] 
for e™ as a function of the true values of |e™ | and 
as well as the true sin 2 2#i3 = 0. This figure 
looks very similar to the corresponding SO-CPV 
figure as a function of sin 2 26*i3 and Sqp '■ There 
is a cutoff at small |e™|, below which the phase 
effects are suppressed, and there is no sensitivity 
close to the CP-conserving solutions 0™ = and 
7r. Therefore, we adopt an approach similar to 
that of SO-CPV. We show in Figure H the frac- 
tion of 0™ (left) and <j>™ (right) for which NSI- 
CPV will be discovered as a function of the |e™ r | 
(left) and |e™ | (right). In this case, the fraction of 
e a/3 represents the stacking of all sensitive regions 
in Figure [1] along any vertical line corresponding 
to any fixed |e™|. In Figure the dependence 
on the true sin 2 26*i3 and <5cp is indicated by the 
shaded regions, whereas the curves correspond to 
the true sin 2 2 9 13 = (i.e., the thick curve in 



4 



the right panel corresponds to the thick curve in 
Figure [J) . 

As it is obvious from the analytical and quan- 
titative discussion in Ref. [25], the v e — > (and 
v e — > v^) appearance channels will dominate the 
determination of 0™ , whereas the — > (and 
— > v^) disappearance channels will dominate 
the the determination of <fi™ T (see also Ref. [33] 
for more analytical discussions). For that reason, 
we obtain a strong dependence on the simulated 
sin 2 2#i3 and Sep for <j>™ (right panel of Figured]), 
because the appearance channels are most sensi- 
tive to these SO parameters, whereas the CP vi- 
olation in </>™ T is hardly affected by these param- 
eters. From Figure (left panel), we can read off 
that NSI-CPV will be discovered for about 80% 
of all possible </>™ T for |e™ | ^0.1 close to the 
current bound. The |e^.| reach is, however, lim- 
ited to |e™ T | > 0.02, which means that any sig- 
nificant improvement of the bound will exclude 
this possibility. For 0™ in the right panel, we 
obtain a picture qualitatively similar to the SO- 
CPV measurement because of the dominance of 
the appearance channels. We obtain a large frac- 
tion of <fi™ of up to 80% in an intermediate range 
0-01 ^ 0-3- From the analytical discus- 

sion in Ref. [25] (Eqs. (10) and (11)), we can- 
not observe the NSI-CPV for too large |e™ |, be- 
cause the terms quadratic in |e™ | will then be too 
large a background for |e™ | > 0(1). This means 
that although the NSI can be easily established, 
the NSI-CPV cannot be established against the 
phase-independent terms. 

4. Mis-interpretation of NSI-CPV 

Let us assume that there is no SO-CPV in na- 
ture, but there is large NSI-CPV. What would 
such a fit look like, and when would one con- 
fuse the NSI-CPV with the SO-CPV? Since 0™ 
is measured in the disappearance channel, which 
is hardly affected by Sep, we focus on in this 
section. 

As an example, we choose a true Sep = (CP 
conservation) and the true </>™ = it/2 (maximal 
CP violation), and we simulate three specific sets 
of true values: 

A sin 2 2<9 13 = 10~ 4 , |e™| = 0.005 



B sin 2 2<9 13 = 10" 4 , |e™| = 0.01 

C sin 2 26»i3 = 1CT 3 , |e™ | = 0.005 

There is some parameter dependence in these 
choices, but the examples are good enough to il- 
lustrate the qualitative main points. We show 
in Figure [31 left three columns, the fits in the 
sin 2 28 13 -5cp (upper row) and |e™|-0™ (lower 
row) planes assuming the NSI scenario includ- 
ing e™. Obviously, in all cases, the minimal 
X 2 = because we simulate the average ex- 
periment performance without statistical fluctu- 
ations. At Ax 2 = 9, NSI-CPV can neither be 
established in case A nor in case C, whereas it 
can be measured in case B (lower row). On the 
other hand, Sqp can be measured in cases B, 
and C, whereas no information can be obtained 
in case A (upper row). Generally speaking, the 
true sin 2 26>i 3 or true |e™ | have to be large enough 
to observe the corresponding phase. In the right 
panel of Figure [3j we illustrate the effect of the 
wrong hypothesis: If only standard oscillations 
are assumed (and |e™ | = in the fit), the min- 
imal x 2 hi all three cases will be non- vanishing, 
and might be confused with statistical fluctua- 
tions. In all three cases, the simulated maximal 
NSI-CPV will be mis-interpreted as SO-CPV, be- 
cause the CP-conserving values Sep — and tt 
can be excluded. Note that especially in cases A 
and C the minimal x 2 is relatively small, and the 
confusion might not be obvious. Even worse, in 
these cases neither SO-CPV nor NSI-CPV would 
be established in the correct fit. 

5. Sensitivity to any CP violation 

Let us suppose that there are NSI, and there 
might be CP violation from either <5cp or c/>™. 
In this case, it may be an interesting question if 
any CP violation can be established, no matter 
of the origin. This question is equivalent to ex- 
cluding the CP conserving solutions (<!>cp,0™) G 
{(0,0),(0,7r),(7r,0),(7r,7r)}. Similar to the "frac- 
tion of Sep" or "fraction of <fi™ T " for the single pa- 
rameter measurement, we define, for a given set 
of true sin 2 2^13 and |e™|, the fraction of (true) 
Sep ® 4^ as the fraction of the Sep d> 4^ T plane 
for which (any) CP violation can be established. 
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No NSI assumed 




Figure 3. NSI-CPV misinterpreted as SO-CPV. The first three columns represent the fits in the sin 2 2#i3- 
<5cp (upper row) and |e™|-(/)™ (lower row) planes including NSI (from e™) for three different sets of true 
values: case A with sin 2 26 13 = 1(T 4 and \e™ \ = 0.005, case B with sin 2 26 13 = 10~ 4 and |e™| = 0.01, 
and case C with sin 2 29i 3 = 10~ 3 and |e™| = 0.005. In all panels, we choose the true Sep = (no 
SO-CPV) and the true 4>™ T = ir/4 (maximal NSI-CPV). In the right panel, we show the same fits in the 
sin 2 26*i3-<5cp plane assuming standard oscillations only (no NSI). In all panels, the contours represent 
A% 2 = 1 (dashed), 4 (solid thin), and 9 (solid thick), and the best-fit values are marked by symbols 
(including the minimum \ 2 )- 



Fraction of Sep <8> for CPV discovery 
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Figure 4. Sensitivity to any (SO or NSI) CPV as 
a function of the true sin 2 26*i3 and true |e™ |. The 
sensitivity is given as a fraction of all (true) Sep 
<j>™ for which CPV can be established (A% 2 = 9). 



This means that for a random pick of Sep and 
<j>™ (uniformly distributed in the phases), this pa- 
rameter will tell the probability that CP violation 
can be established. We show this sensitivity to 
any (SO or NSI) CPV as a function of the true 
sin 2 26*13 and true |e™| in Figure |U Obviously, if 
both sin 2 26 13 and |e™ | are small, no CPV can be 
established. If only one of the parameters is small, 
the sensitivity corresponds to the single parame- 
ter measurement (for the other parameter and its 
phase being marginalized over). If both sin 2 26>i3 
and |e™ | are large and in their optimal ranges, the 
fraction of <5cp <8> 4>™ can reach more than 95%. 
In this case, CPV cannot be established in only 
very small regions around the four CP-conserving 
solutions given above. In the probability inter- 
pretation, the chance to find CPV is therefore, 
in fact, larger if both large sin 2 26*i 3 and |e™ | are 
present, compared to the case of no NSI. 
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6. Summary and conclusions 

We have demonstrated that, in a neutrino fac- 
tory corresponding to the current IDS-NF base- 
line setup, CP violation from non-standard neu- 
trino interactions in matter could be measured if 
|e™ I > 0.02 or |e™| > 0.005. This observation is 
almost independent of the true values of sin 2 2$i3 
and <5cp. We have also shown that there are re- 
gions in parameter space where an erroneous as- 
sumption of standard oscillations only will lead to 
a fake CP violation signal in the presence of NSI, 
even if CP is conserved in Sep- Finally, we have 
discussed the chances to measure any CP viola- 
tion in the presence of both large enough sin 2 2#i3 
and |e™ |. We have found that CP violation will 
be found in up to 95% of all possible phase combi- 
nations for Sep and 0™. For the quantification, 
we have used performance indicators similar to 
the standard CP violation measurement. 

We conclude that even if there is no CP vio- 
lation in standard oscillations or sin 2 26i% is too 
small to detect it, a neutrino factory has the 
chance to find CP violation from new physics ef- 
fects. In any future analysis, it is therefore impor- 
tant to carefully consider the possibility of non- 
standard effects, in order not to overlook an even 
more interesting hint for new physics. 
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